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Q , Abstract 

O , 

| Spectral efficiency for asynchronous code division multiple access (CDMA) with random spreading is calculated 

in the large system limit allowing for arbitrary chip waveforms and frequency-fiat fading. Signal to interference and 

O ! 

noise ratios (SINRs) for suboptimal receivers, such as the linear minimum mean square error (MMSE) detectors, are 
derived. The approach is general and optionally allows even for statistics obtained by under-sampling the received 



00 _ 
£N| | signal. 

All performance measures are given as a function of the chip waveform and the delay distribution of the users in 
the large system limit. It turns out that synchronizing users on a chip level impairs performance for all chip waveforms 
c/5 \ with bandwidth greater than the Nyquist bandwidth, e.g., positive roll-off factors. For example, with the pulse shaping 
demanded in the UMTS standard, user synchronization reduces spectral efficiency up to 12% at 10 dB normalized 
signal-to-noise ratio. The benefits of asynchronism stem from the finding that the excess bandwidth of chip waveforms 
k>" ' actually spans additional dimensions in signal space, if the users are de-synchronized on the chip-level. 
qq | The analysis of linear MMSE detectors shows that the limiting interference effects can be decoupled both in the 

user domain and in the frequency domain such that the concept of the effective interference spectral density arises. 
This generalizes and refines Tse and Hanly's concept of effective interference. 

In Part II, the analysis is extended to any linear detector that admits a representation as multistage detector and 
guidelines for the design of low complexity multistage detectors with universal weights are provided. 
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I. Introduction 

The fundamental limits of synchronous code-division multiple-access (CDMA) systems and the loss in- 
curred by the imposition of suboptimal receiving structures have been thoroughly studied in different scenar- 
ios and from different perspectives. On the one hand, significant efforts have been devoted to characterize 
the optimal spreading sequences and the corresponding capacities 0J, 0, Q. On the other hand, very in- 
sightful analysis 01, 0, (6l, (7]|, jHl resulted from modelling the spreading sequences by random sequences 
flU. In fact, as both the number of users K and the spreading factor N tend to infinity with a fixed ratio, 
CDMA systems with random spreading show self-averaging properties. These enable the description of the 
system in terms of few macroscopic system parameters and thus provide a deep understanding of the system 
behavior. 

In the literature, the fundamental limits of CDMA systems and the asymptotic analysis of linear multiuser 
detectors under the assumption of random spreading sequences is overwhelmingly focused on synchronous 
CDMA systems. While the assumption of user synchronization allowed for accurate large-system analysis, 
it is not realistic for the received signal on the uplink of a cellular CDMA system, in particular if users 
move and cause varying delays. Therefore, it is of theoretical and practical interest to extend the analysis of 
CDMA systems with random spreading to asynchronous users. This holds in particular, as we will see that 
asynchronous users are beneficial from a viewpoint of system performance. 

The analysis of asynchronous CDMA systems using a single-user matched filter as receiver was first given 
in [[Toll . A rich field of analysis of asynchronous CDMA systems with conventional detection at the receiver 
is based on Gaussian approximation methods. An exhaustive overview of these approaches exceeds the 
scope of this work, which is focused on the analysis of asynchronous CDMA systems with optimal joint 
decoding or linear multiuser detection. The interested reader is referred to 0TJ |[T2| and references therein 
for asynchronous CDMA with single-user receivers. 

The analysis and design of asynchronous CDMA systems with linear detectors is predominantly restricted 
to consider symbol- asynchronous but chip-synchronous signals, i.e., the time delays of the signals are mul- 
tiples of the chip interval. The effect of chip-asynchronism is eventually analyzed independently iTOl . In 
this stream are works that optimize the spreading sequences to maximize the sum capacity lTl4ll and analyze 
the performance of linear multiuser detectors flS, HJl, OH, EL In US, OS, the linear MMSE detec- 
tor for symbol-asynchronous but chip-synchronous systems is shown to attain the performance of the linear 
MMSE detector for synchronous systems as the size of the observation window tends to infinity by empirical 
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and analytical means, respectively. However, they verify numerically that the performance of linear MMSE 
detectors is severely impaired by the use of short observation windows. Additionally, lfl"5l provides the 
large-system SINR for a symbol whose chips are completely received in an observation window of length 
equal to the symbol interval T s . In IfToll . IfTTl . the analysis of linear multistage and MMSE detectors is ex- 
tended to observation windows of arbitrary length. Furthermore, a multistage detector structure that does 
not suffer from windowing effects and performs as well as the multistage detector for synchronous systems 
is proposed. 

In |fT3ll . [fT8l . the effects of chip asynchronism are analyzed assuming bandlimited chip pulses. In [18], the 
chip waveform is assumed to be an ideal Nyquist sine function, i.e., a sine function with bandwidth equal 
to half of the chip rate. The received signal is filtered by a lowpass filter (or, equivalently, a filter matched 
to the chip waveform) and subsequently sampled at the time delay of the signal of the user of interest with 
a frequency equal to the chip ratdj. Reference |[T8l proves that the SINR at the output of the linear MMSE 
detector converges in the mean-square sense to the SINR in an equivalent synchronous system. In lfl"3l the 
wider class of chip pulses which are inter-chip interference free at the output of the chip matched filter is 
considered. In the following we will refer to this class of chip pulses as square root Nyquist chip pulses. 

In |[T9ll , ll20ll . the performance of the linear MMSE detector with completely asynchronous users and 
chip waveforms limited to a chip interval is analyzed. However, the observation window in lfT9l spanned 
only a single symbol interval not yielding sufficient discrete-time statistics; the resulting degradation in 
performance was pointed out later in lfT3l . lfl"5l . 

As discussed above, previous approaches to the analysis of asynchronous CDMA with multiuser detection 
were only concerned with, if and how asynchronism can be prevented from causing performance degrada- 
tion. However, asynchronism is known to be beneficial for CDMA systems with demodulation by single-user 
matched filters (e.g., [10]). One of the main contributions of this paper is to show that benefits from asyn- 
chronism are not inherent to single-user matched filters but they are a general property of CDMA systems. 
We quantify those benefits in terms of spectral efficiency and SINRs in the large-system limit. 

Compared to synchronous systems, the analysis of asynchronous CDMA raises two additional issues: 
(i) the way statistics are formed, trading complexity against performance, and (ii) the effects of excess 
bandwidth, chip-pulse shaping and the users' delay distribution. 

The optimum multiuser detector in ll2lTl is based on the sufficient statistics obtained as output samples of 

1 The chip rate satisfies the condition of the sampling theorem in this case. 
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a bank of filters matched to the symbol spreading waveforms of all users. The decorrelating detector in ll22ll 
and the linear MMSE detector in 11231 benefit from the same sufficient statistics. A method to determine 
the eigenvalue moments of a correlation matrix in asynchronous CDMA systems using such statistics and 
square root Nyquist chip pulse waveforms is proposed in Il24ll . 

An alternative approach to generate useful statistics, which in general are not sufficient, is borrowed from 
synchronous systems. The received signals are processed by a filter matched to the chip waveform and 
sampled at the chip rate. This approach is optimum for single-user communications and chip-synchronous 
multi-user communications, but causes aliasing to the signals of de-synchronized users if the chip wave- 
form has non-zero excess bandwidth. Discretization schemes using chip matched filters and sampling at 
the Nyquist rate are studied in G51 . There, the notion of approximate sufficient statistics was introduced. 
Furthermore, conventional CDMA systems with chip waveforms that approximate sine pulses were shown 
to outperform systems using rectangular pulse shaping. Moreover, it was conjectured that sine pulses are 
optimal for CDMA systems with linear MMSE multiuser detection. 

In systems with bandlimited waveforms, sampling at a rate faster than the Nyquist rate leads to the same 
performance as the optimal time-discretization proposed in |[2T|. |[22|. Il23l if the condition of the sampling 
theorem is satisfied lfl"3~1 . In contrast to the bank of symbol matched filters in ETI . this approach has the 
advantage that the time delays of the users' signals need not be known before sampling. 

The impact of the shape and excess bandwidth of the chip pulses received attention in ll26l . Il27l . G51 . 
In 11261 . Il27l an algorithm for the design of chip-pulse waveforms for CDMA systems with conventional 
detection has been proposed. The design criterion consists of minimizing the bit error rate at the output of 
a single user matched filter in asynchronous CDMA systems while enforcing certain constraints on the chip 
waveforms. 

This work is organized in seven additional sections. Section UH gives a brief overview of the main re- 
sults found in this work. Sections [TIT] and [IV] introduce notation and the system model for asynchronous 
CDMA, respectively. Section [V] focuses on the analysis of linear MMSE detectors and introduces the main 
mathematical tools for analysis of the fundamental limits of asynchronous CDMA. In Section [VH the spec- 
tral efficiency of optimal joint decoding is derived on the basis of the results for the linear MMSE detector 
exploiting the duality between mutual information and MMSE. Section [VTT1 addresses the extension of the 
presented results to more general settings. Some conclusions are drawn in Section IVIIII 
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II. Main Results 

Before going into the main results of this work, it is helpful to get some intuition on asynchronous CDMA 
systems. First of all, one might be interested in the question which chip waveform gives the highest spectral 
efficiency for otherwise arbitrary system parameters, like pulse shape, pulse width in time and frequency, 
system load, etc. There is a surprisingly easy answer to this question that does not require any sophisticated 
mathematical tools: 

Proposition 1 Without constraints on the system parameters, Nyquist sinc-pulses maximize spectral effi- 
ciency . 

Proof: The proof is by contradiction. First, it is well-known that the spectral efficiency of a single user 
channel is maximized by Nyquist sinc-pulses. Further, we know from [21 that the spectral efficiency of a 
synchronous CDMA system with Nyquist sinc-pulses becomes identical to the spectral efficiency of a single 
user channel, as the load converges to infinity. Finally, the multiuser system can never outperform the single 
user system, since we could otherwise improve a single user system by virtually splitting the single user into 
many virtual users. Thus, the Nyquist sine pulse is optimum also for the multi-user system. ■ 

Note that, from the previous proof, the Nyquist sine pulse is optimum for an infinite system load. However, 

we cannot judge whether the optimum is unique from the line of thought proposed in our proof. In fact, a 

straightforward application of a more general result in this paper (shown in the Appendix© is the following: 



Proposition 2 Asynchronous CDMA systems with any sinc-pulses, no matter whether they are constrained 
to the Nyquist bandwidth or to a larger, or even to a smaller bandwidth^ and users whose empirical delays 
are uniformly distributed within a symbol interval achieve the same spectral efficiency as a single user 
channel, if the load converges to infinity. 

The optimization of the system load neither gives the theoretically most interesting cases to consider 
nor the practically most relevant. Let us, thus, look at which chip waveforms achieve the highest spectral 

2 It is worth noticing that any chip pulse with bandwidth smaller than the Nyquist bandwidth necessarily introduces inter-chip interference. 
November 28, 2009 
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Fig. 1 

Decomposition of a spreading sequences of length N = 12 into M = 3 virtual spreading sequences. 

efficiency for a fixed load. Surprisingly, the result is not a little bit more useful for practical applications: 

Corollary 1 The chip waveform that maximizes the spectral efficiency for a given finite load and given chip 
rate has vanishing bandwidth. Furthermore, the maximum spectral efficiency is the same as the one for the 
single user channel. 

Proof: The corollary follows directly from Proposition^ Note that Proposition[2]holds for an arbitrary 
chip rate and an arbitrary bandwidth of the chip waveform and states that the single user bound is reached 
at infinite load. Though, the corollary is stated for a given finite load, we are free to decompose each 
physical user into M virtual users and let M increase to infinity such that the virtual load becomes infinite. 
Therefore, we take a user's signal and divide it into M data streams that are time-multiplexed in such a 
way as to result in the same physical transmit signal for that user. An example of such a decomposition 
is illustrated in Fig. [Q Applying this idea to each of the K physical users, we have created MK virtual 
users. Furthermore, the chip interval has grown from T c to MT C and the virtual users are asynchronous 
with a discrete uniform distribution of delays within the virtual chip interval of length MT C . Consider now 
a sine pulse of bandwidth 1 / (2MT C ) as chip waveform. If we take the limit M — > oo for the system of 
virtual users, the delay distribution converges to the uniform distribution within the virtual chip interval and 
the number of virtual users converges to infinity. Thus, Proposition [2] applies and the single user bound is 
reached. Therefore, this choice of chip waveforms whose bandwidth vanishes is optimal. ■ 
Optimizing chip waveforms to maximize spectral efficiency has proven to hardly aid the practical design 
of CDMA systems, since the optima are achieved for system parameters, e.g., infinite load and/or vanish- 
ing bandwidth, that are far from the limits of practical implementation. Furthermore, the choice of the 
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chip waveform is influenced by many other factors than spectral efficiency like the difficulty to implement 
steeply decaying frequency filters and the need to keep the peak-to-average power ratio of the continuous- 
time transmit signal moderate. Therefore, many commercial CDMA systems, e.g., the Universal Mobile 
Telecommunication System (UMTS), use chip waveforms with excess bandwidth. The UMTS standard 
uses root-raised cosine pulses with roll-off factor 0.22. Motivated by the theoretical findings above and the 
practical constraints on the design of chip waveforms, the rest of this paper puts the focus on the perfor- 
mance analysis of CDMA system with a given fixed chip waveform. As it will be seen, this gives rise to a 
rich collection of insights into CDMA systems with asynchronous users. The main results are summarized 
in the following. 

CDMA systems using chip pulse waveforms with bandwidth B not greater than half of the chip rate Jr, 

1 c 

i.e., B < , perform identically irrespective of whether the users are synchronized or not for a large class of 
performance measures. Furthermore, our result generalizes the equivalence result for the ideal Nyquist sine 
waveform (with bandwidth ^r) in |fT8l to any chip pulse satisfying the mentioned bandwidth constraint and 
to any linear multistage detector and the optimal capacity-achieving joint decoder. Note that the performance 
is independent of the time delay distribution. Increasing the bandwidth of the chip waveform above ^r, i.e., 
allowing for some excess bandwidth as it is customary in all implemented systems, the behaviors of CDMA 
systems change substantially. They depend on the time delay distribution and the equivalence between 
synchronous and asynchronous systems is lost. 

For any choice of chip waveform, we capture the performance of a large CDMA system with linear 
MMSE detection by a positive definite frequency-dependent Hermitian matrix Y(f2) whose size is the ratio 
of sampling rate to chip rate (the sampling rate is a multiple of the chip rate). We require neither the absence 
of inter-chip interference, nor that the samples provide sufficient statistics, nor a certain delay distribution. 
Unlike for synchronous users, the multiuser efficiency |[2~8l in the large-system limit is not necessarily unique 
for all users. The matrix Y(f2) reduces to a scalar frequency dependent function t]{oj) in cases where 
oversampling is not needed. Interestingly, the same holds true even in cases with excess bandwidth if the 
delay distribution is uniform. The scalar r/(a>) can be understood as a multiuser efficiency spectral density 
with the multiuser efficiency being its integral over frequency oj. We find that in large systems, the effects 
of interference from different users and interference at different frequencies decouple. We, thus, generalize 
Tse and Hanly's flU concept of effective interference to the concept of effective interference spectral density 
which decouples the effects of interference in both user and frequency domain. 
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Excess bandwidth can be utilized if users are asynchronous. While excess bandwidth is useless for syn- 
chronized systems in terms of multiuser efficiency, i.e., all square root Nyquist pulses perform the same 
regardless of their bandwidth, desynchronizing users improves the performance of any system with non- 
vanishing excess bandwidth. 



III. Notation and Some Useful Definitions 



Throughout this work, upper and lower boldface symbols are respectively used for matrices and vectors 
spanning a single symbol interval. Matrices and vectors describing signals spanning more than a symbol 
interval are denoted by upper boldface calligraphic letters. 

In the following, we utilize unitary Fourier transforms both in the continuous time and in the discrete time 
domain. The unitary Fourier transform of a signal s(t) in the continuous time domain is given by S(u) = 
-|= s(t)e~ JW *dt. The unitary Fourier transform of a sequence {. . . , c_i, c , ci, . . .} in the discrete time 
domain is given by c(Q) = -j= ^^=-00 c n^~^ n - We will refer to them shortly as Fourier transform. 
Throughout this work, u and Q denote the angular frequency and the angular frequency normalized to the 
chip rate, respectively. A function in Q has support in the interval [— tt, n], or translations of it. 

For further studies it is convenient to define the concept of r-block-wise circulant matrices of order N: 



Definition 1 Let r and N be positive integers. An r-block-wise circulant matrix of order N is an rN x N 
matrix of the form 



( 



C 



B Q B 1 B N -i 
Bn-i B ■ ■ ■ Bn-2 



\ 



\ B 1 B 2 



(1) 



B J 



with Bi = (c M , c 2 ,i, ...,Cr 



In the matrix C, an r x N block row is obtained by a circular right shift of the previous block. Since the 
matrix C is univocally defined by the unitary Fourier transforms of the sequences 



JV-l 



c s {n) 



k=0 



Cs ,k 



-jUk 



r 

3***35 



(2) 
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there exists abijection $ from the frequency dependent vector c(Q) = [ci(f2), C2(fi), . . . , c r (f2)] to C. Thus, 

d{c(n)} = c. O) 

Furthermore, the superscripts - T and - H denote the transpose and the conjugate transpose of the matrix 
argument, respectively. I n is the identity matrix of size n x n and C, Z, Z + , N, and R are the fields of 
complex, integer, nonnegative integer, positive integer, and real numbers, respectively, tr(-), || ■ || , and | • | are 



the trace, the Frobenius norm, and the spectral norm of the argument, respectively, i.e., || A|| = y tr(AA^), 
|A| = max x^AA^x. diag(-) : C n i— > C nxn transforms an n-dimensional vector into a diagonal matrix of 

X H X<1 

size n x n having as diagonal elements the components of the vector in the same order. E{-} and Pr{-} are 
the expectation and probability operators, respectively. 5ij is the Kronecker symbol and 5(A) is Dirac's delta 
function. X = (£yOi=iV."ni * s tne ni x 77-2 matr i x whose (i, j)-element is the scalar x^-. X = (-X^j)i=i'"'ni 
is the ni^i x n 2 q2 block matrix whose (z,^) -block is the qi x q 2 matrix X^. The notation [-J is adopted for 
the operator that yields the maximum integer not greater than its argument and xmody denotes the modulus, 



i.e., xmody = x 



y. Furthermore, x( x -A) denotes the indicator function of the variable x on the 



set A and x( x £ A) = lif x G A and zero otherwise. 



IV. System Model 

Let us consider an asynchronous CDMA system with K users in the uplink channel. Each user and the 
base station are equipped with a single antenna. The channel is flaj^] fading and impaired by additive white 
Gaussian noise. Then, the signal received at the base station, in complex base-band notation, is given by 

K 

y{t) = ^2a k s k (t - T k ) + w(t) tG(-oo,+oo). (4) 

k=l 

Here, a k is the received signal amplitude of user k, which takes into account the transmitted amplitude, the 
effects of the flat fading, and the carrier phase offset; r k is the time delay of user k; w(t) is a zero mean 
white, complex Gaussian process with power spectral density N ; and s k (t) is the spread signal of user k. 
We have 

+oo 

»*(*)= E MH4 m) (*)> (5) 

m=—oo 

3 Flat fading is no restriction of generality here as long as the excess delay is much smaller than the symbol interval T B . This is, as the effect 
of multi-path can be incorporated into the shape of the chip waveform. 
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where bk[m] is the m th transmitted symbol of user k and 



iV-l 



4 (*) = s k,m[n}ip(t - mT s - nT c ) 



(6) 



n=0 



is its spreading waveform at time m. Here, S ™ is the spreading sequence vector of user k in the m 
symbol interval with elements Sk, m [n], n — 0, . . . , N — 1. T s and T c = ^ are the symbol and chip interval, 
respectively. 

The users' symbols bk\m] are independent and identically distributed (i.i.d.) random variables with 
E{|6fc[m]| 2 } = 1 and E{&fc[m]} = 0. The elements of the spreading sequences 4 N ^ assumed to 
be i.i.d. random variables with E{|sfc m [?T,]| 2 } = i and E{s/- jm [n]} = 0. This assumption properly models 
the spreading sequences of some CDMA systems currently in use, such as the long spreading codes of the 
FDD (Frequency Division Duplex) mode in the UMTS uplink channel. 

The chip waveform tf)(t) is limited to bandwidth B and energy = J_°° |^(t)| 2 dt. Because of the 
constraint on the variance of the chips, i.e., E{|s fcm [/c]| 2 } = the mean energy of the signature waveform 



1) user 1 as reference user so that T\ = 0, 

2) the users are ordered according to increasing time delay with respect to the reference user, 

3) the time delay to be, at most, one chip interval so that r fc G [0, T c ). 

Assumptions 1 and 2 are without loss of generality ll2~9l . Assumption 3 is made for the sake of clarity and it 
will be removed in Section [VlT] where the results are extended to the general case with G [0, T s ). 

At the receiver front-end, the base band signal is passed through a filter with impulse response g(t) and 
corresponding transfer function G{uo) normalized such that f*™ \g(t)\ 2 dt = 1. We denote by (f>(t) the 
response of the filter to the input i/i(t), i.e., 4>{t) = g(t) * tp(t) and by its Fourier transform. The filter 
output is sampled at rate with r G N. For further convenience, we also define = \<p(t)\ 2 dt. 

Throughout this work we assume that the filtered chip pulse waveform cf){t) is much shorter than the sym- 
bol waveform, i.e., 4>{t) becomes negligible for \t\ > t and t T s . This technical assumption is usually 
verified in the systems with large spreading factor we are considering. It allows to neglect intersymbol inter- 
ference. Thus, focusing on a given symbol interval, we can omit the symbol index m and the discrete-time 
signal at the front-end output is given by 




K 




(7) 



k=i 
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with sampling time p G {0, . . . , rN— 1} and 



N-l 



c k(t) = ^2 S k[ n ]<t> (* ~~ nTc ) ■ 



(8) 



n=0 



Here, w[p] is discrete-time, complex- valued noise. In general, w[p] is not white, although the continuous 
process was white. However, it is white, if g(t) * g(—t) is Nyquist with respect to the sampling rate. In this 
latter case, the noise variance is a 2 = This expression accounts for the normalization of the front end 

J- c 

filter. 

In order to cope with the effects of oversampling, we consider an extended signal space with virtual 
spreading sequences of length rN. The virtual spreading sequence of user k is given by the iVr-dimensional 
vector 

$ k s k (9) 



where si 



f Mo] 



s k [N-l]) T 
( 



^r k 



0- 7 



-Tk 



■ ■ 4>(-N+2)T c -r k 



(10) 



is a Nr x N block matrix taking into account the effects of delay and pulse shaping. Its blocks are the 
vectors <p x = (J>(x), (p (x + ^) , . . . , (x + r -^-T c )) T . In that way, we have described user k's, continuous- 
time channel with continuous delays canonically by the discrete-time channel matrix Note that 
solely depends on the delay of user k, the oversampling factor r, the chip waveform, and the receive filter. 

Structuring the matrix in blocks of dimensions r x 1, it is block- wise Toeplitz. As well known ll30l . 
OTI . block- Toeplitz and block-circulant matrices are asymptotically equivalent in terms of spectral distribu- 
tion. This asymptotic equivalence is sufficient for us, since in this work we focus on performance measures 
of CDMA systems which depend only on the asymptotic eigenvalue distribution. Similar asymptotically 
tight approximations are used in the large system analysis of CDMA in frequency-selective fading ll32l . 
031,1133. 

The equivalent block-circulant matrix is given by 



(r-l)T c 
r 



(ID 
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where 

, +00 

m r) = - £ e^ (n+2 ^$* + 27TI/)) (12) 



1/=— 00 



is the spectrum of the chip waveform 0(t) delayed by r and sampled at rate 1/T C . Thus, we replace the 
block- Toeplitz matrix for our asymptotic analysis by the block-circulant matrix <frfc in the following and 
use the virtual spreading sequences 

v k = ® k s k . (13) 

Let S be the rN x K matrix of virtual spreading, i.e., S = (3>iSi, $2^2, • • • &k&k), A the. K x K 
diagonal matrix of received amplitudes, H = SA, and b and 

y = Hb + w (14) 

the vectors of transmitted and received signals, respectively. Additionally, h k denotes the k th column of the 
matrix H. Finally, we define the correlation matrices T = HH H , R = H H H and the system load (3 = ^. 

V. Linear MMSE Detection 

The linear MMSE detector d k generates a soft decision b k = d^y of the transmitted symbol b k based on 
the observation y. It can be derived from the Wiener-Hopf theorem [[331 and is given by 

d k = V{yy H }- l V{bly} (15) 

with the expectation taken over the transmitted symbols b and the noise. Specializing the Wiener-Hopf 
equation to the system model CEH) yields 

d k = (HH H + a 2 I)~ l h k (16) 
= c-{H k H% + a 2 I)- l h k (17) 

for some cel. Here, H k is the matrix obtained from H suppressing column h k . The second step follows 
from the matrix inversion lemma. 

The performance of the linear MMSE detector is measured by the signal-to-interference-and-noise ratio 
at its output Il28l 

SINR fc = h% (H k H% + a 2 I)~ l h k . (18) 
The SINR can be conveniently expressed in terms of the multiuser efficiency r] k [|28l 

SINR fc = Mr%*- (19) 
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The multiuser efficiency is a useful measure, since for large systems it is identical for all users in special 
cases and it is related to the spectral efficiency 0, ll36ll . Il37i 

The SINR depends on the spreading sequences, the received powers of all users, the chip pulse shaping, 
and the time delays of all users. To get deeper insight on the linear MMSE detector it is convenient to analyze 
the performance for random spreading sequences in the large system limits, i.e., as K, N — > oo with the ratio 
jf — > (3 kept fixed. The large-system analysis will identify the macroscopic parameters that characterize a 
chip-asynchronous CDMA system and the influences of chip pulse shaping and delay distribution. 

In this section we present the large system analysis of a linear MMSE detector for chip-asynchronous 
CDMA systems with random spreading. Provided that the noise at the output of the front end is white, 
the analysis applies to CDMA systems using either optimum or suboptimum statistics, any chip pulse wave- 
forms, and any set of time delays in [0, T c ) if their empirical distribution function converges to a deterministic 
limit. 

In Appendix UH we derive the following theorem on the large-system performance of chip-asynchronous 
CDMA: 

Theorem 1 Let A G <C KxK be a diagonal matrix with k th diagonal element a k G C and T c a positive 
real. Given a function $(u;) : R — > C, let <p(Q, r) be as in 4721) . Given a positive integer r, let <fr fc , 
k = 1, . . . , K, be r -block-wise circulant matrices of order N defined in 4771) . Let H = SA with S = 
[<&iSi, 3> 2 s 2 , • • • , $ks k ] and s k G C^ 1 . 

Assume that the function \ $ (u) \ is upper bounded and has finite support. The receive filter is such that the 
sampled discrete-time noise process is white. The vectors s& are independent with i.i.d. circularly symmetric 
Gaussian elements. Furthermore, the elements a k of the matrix A are uniformly bounded for any K. The 
sequence of the empirical joint distributions F^l T (A,r) = ^^2 k =iX(^ > \ a k\ 2 )x( T > r k ) converges 
almost surely, as K — > oo, to a non-random distribution function F\a\ 2 ,t(\ t )- 

Then, given the received power laJ 2 , the time delay r k , and the variance of the white noise a 2 = 
the SINR of user k at the output of a linear MMSE detector for a CDMA system with transfer matrix H 
converges in probability as K, N — >• oo with ^ — > (3 and r fixed to 
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lim SINR fc 

K=/3N-^oo 



|QfcT 

2vr 



Aj r (ft,T fe )Y(ft)A 0ir .(ft,Tfc)dft 



(20) 



where Y(ft) is the unique positive definite r x r matrix solution of the fixed point matrix equation^ 

— 7T < ft < 71, 



T _, = ^ + /* A^(»,,)Ag r («,r)^ a , T (A,r) 

•/ l + i/: 7r Aj r (fi',r)T(n')A 0)r (O',r)dn' 



(21) 



and 



/ 



A^ r (fi,r) 



0(ft,r 
0(ft,r- 



T c (r-1)- 



(22) 



The performance of the linear MMSE detector operating on not necessarily sufficient statistics is com- 
pletely characterized by 

1) an r x r matrix-valued transfer function Y(ft) and 

2) the frequency and delay dependent vector A<£ iT .(ft, r) 

The multiuser efficiency varies from user to user and depends on the time delay of the user of interest only 
through A^, r (ft, Tfe). We can define an SINR spectrum 



SINR fc (ft) 



|QfcT 
2tt 



Ag r (fi,7*)T(fi)A*, r (fl,7*) 



(23) 



in the normalized frequency domain — 7r < ft < 7r, or, equivalently, a spectrum of the multiuser efficiency. 
The system performance is in both cases obtained by integration over the spectral components. 

The fixed point equation (121")) clearly reveals how and why synchronous users are the worst case for a given 
chip waveform. We know from 11371 that to each large multiuser system, there is an equivalent single user 
system with enhanced noise, but otherwise identical performance. In the present case with oversampling 
factor r, the equivalent single user system is a frequency- selective MIMO (multiple-input multiple-output) 
system with r transmit and r receive antenna and governed by the r x r channel transfer matrix 

AA^ r (ft, r) Aj r (ft, r)d*j A | a>T (A, r) 



P ' l + A^Aj r (ft',r)T(ftOA,, r (ft',r)dft'- 

4 Here, the integration measure is meant to denote dFj^p^A, r) = /|a| 2 ,t(A, i")dAdr in case such a representation exists. 
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Note that this matrix is an integral of an outer product over the delay distribution. Thus, for constant delay, 
i.e., chip-synchronization, the matrix has rank one. No additional dimensions in signal space can be spanned. 
For distributed delays, the rank of the matrix can be as large as the oversampling factor r. Driving the 
equivalence even further, the equivalent MIMO system can be transformed into an equivalent CDMA system 
with spreading factor r and spreading sequences A^ r (f2, r^). In this model, equal delays in the real CDMA 
system correspond to users with identical signature sequences in the equivalent CDMA system. 

One cannot increase performance unboundedly by faster oversampling, as not all modes of the equivalent 
r-dimensional MIMO system can be excited with a chip waveform of limited excess bandwidth due to the 
projection onto the spectral support of the chip waveform in (|2Q|) . In order to utilize the excess bandwidth 
of the system, we need two ingredients: 1) Time delays separating the users by making the signatures in 
the equivalent system differ. 2) A receiver that transforms the continuous-time receive signal into sufficient 
discrete-time statistics, e.g., by oversampling. A lack of different delays leads to a system where only a 
single eigenmode of the equivalent MIMO system is excited. A lack of oversampling leads to a system 
where more eigenmodes are excited, but are not converted into discrete time. 

Additional intuitive insight into the behavior of the asynchronus CDMA systems can be gained by focus- 
ing on CDMA systems with uniformly distributed delay. In this case, Theorem \T\ can be formulated with a 
single scalar fixed point equation by moving from the frequency Cl that is normalized to the chip rate to the 
unnormalized frequency oj. This yields the following corollary: 

Corollary 2 Let us adopt the same definitions as in Theorem [7J and let the assumptions of Theorem [7J be 
satisfied. Additionally, assume that the random variables A and r in F\^2^{\ r) are statistically indepen- 
dent and the random variable r is uniformly distributed in [0, T c ). Furthermore, let $(o;) vanish outside 
the interval [—2nB; +2irB] with B < ^r. Then, the multiuser efficiency of the linear MMSE detector for 
CDMA converges in probability as K, N — > oo with jt — > (3 and r fixed to 

+2ttB 
1 f 

lim r]k = r] = — / 77 (u) duj (25) 



K=f3N^oo 2tT 

-2-kB 



where the multiuser efficiency spectral density rj (u) is the unique solution to the fixed point equation 

1 (3 f\dF ]A] 2{X) 



(26) 
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and is zero for \oj\> 2nB. 

Theorem [T]is specialized to Corollary [2] in Appendix ITlTl 

Under the conditions of Corollary [2] the multiuser efficiency of the linear MMSE detector in asynchronous 
systems is the same for all users. 

Rewriting (1251) and (T2~6l) in terms of SINRs, these equations can be interpreted similarly to the correspond- 
ing equations in [[51 for synchronous systems when the concept of effective interference is generalized to the 
concept of effective interference spectral density. Let P(u, A) = ^- 1$ (u;)| 2 be the power spectral density of 
the received signal for a user having received power A. Then, the result in Corollary [2] can be expressed as 

SINR fc = — f smr fc (^)dcu (27) 



2tt 

-2-kB 

where the SINR spectral density sinrfc(u;) is given by 

• / x = Pi^ W\ 2 ) r2R , 

k{ ) Nq + pEx{I{p{uAakl 2^p {ujA)ySmRk)} ^ > 

with the effective interference spectral density 

I(P(u, \a k n P( W> A), SINR fc ) = ^flf^Lr, • (29) 

P{u, \a k \ z ) + r (a;, ApUNRfc 

Heuristically, this means that for large systems the SINR spectral density is deterministic and given by 

smr fc (a;) w Ar ■ , — — — — . — , . — — . — , . nTATT ^ . . (30) 







2 ) 




a k \ 2 ),P(u, 




2 ),SINR fc ) 



This result yields an interpretation of the effects of each of the interfering users on the SINR of user k 
similar to the case of synchronous systems in Q. The impairment at frequency u can be decoupled into 
a sum of the background noise and an interference term from each of the users at the same frequency. 
The cumulated interference spectral density at frequency u depends only on the received power density 
of the user of interest at this frequency, the received power spectral density of the interfering users at this 
frequency, and the attained SINR of user k. In other words, in asynchronous systems we have a decoupling 
of the effects of interferers like in synchronous systems and an additional decoupling in frequency. The 
term I(P(u, \a k \ 2 ), P(uj, |a., | 2 ), SINR fc ) is the effective interference spectral density of user j onto user k at 
frequency uj for a given SINR of user k. 

Sine waveforms have a particular theoretical interest. In the following we specialize Corollary [2] to this 
case. 
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Corollary 3 Let us adopt the definitions in Theorem [7] and let the assumptions of Corollary [2] be satisfied. 
Given a positive real a, we assume that 



't,- : 



for lfl<^, 

(3D 

otherwise. 



corresponding to a sine waveform with bandwidth B = ^fr and unit energy. Then, the multiuser efficiency 
of the linear MMSE detector converges in probability as K, N — > oo with — > (3 to 

lim T] k = ?7 s i nc (32) 

K=f3N^oo 

where the multiuser efficiency r] sinc is the unique positive solution to the fixed point equation 

1 f AdF.WA) 

1 + - / Ar 1 v ; . (33) 



/ N + \r] 

s 

We recall that the multiuser efficiency of a linear MMSE detector for a synchronous CDMA system 



satisfies llSl 



1 f \dF\*a(\) 

1 + /3 AT 1 V ; . (34) 



Vsyn J N Q + \T) syQ 

This result holds for synchronous CDMA systems using any chip pulse waveform with bandwidth B > -^r 
and satisfying the Nyquist criterion. Thus, it also applies to sine pulses whose bandwidth is an integer 
multiple of Then, Corollary [3] shows the interesting effect that an asynchronous CDMA system using a 
sine function with bandwidth B = ^fr as chip pulse waveform performs as well as a synchronous CDMA 
system with bandwidth ^r, r G N, and system load 0' = ^. This implies that only asynchronous CDMA 
has the capability to trade the excess bandwidth of the chip pulse waveform against the spreading factor 
while synchronous CDMA has not. In other words, asynchronous CDMA offers to trade degrees of freedom 
in the frequency domain provided by the excess bandwidth of the chip pulse waveform against degrees of 
freedom in the time domain provided by spreading. 

This phenomenon is similar to the resource pooling in CDMA systems with spatial diversity discovered 
in ll38l . There, the degrees of freedom in space provided by multiple antennas at the receiver could be traded 
against degrees of freedom in time provided by the spreading. In order to make resource pooling happen, it 
is necessary that the steering vectors of the antenna arrays point into different directions. This condition is 
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equivalent to requiring de-synchronization among users. If all users experience the same delay, this is like 
having totally correlated antenna elements. 

In Corollary [31 the bandwidth of the sine waveform may be either larger or smaller than the Nyquist band- 
width. For larger bandwidth, we get a resource pooling effect, for smaller bandwidth we create inter-chip 
interference and what could be called anti-resource pooling. Inter-chip interference is no particular cause of 
concern. In contrast, the effect of anti-resource pooling is to virtually increase the load, i.e., squeezing the 
same number of data into a smaller spectrum is equivalent to squeezing more users into the same spectrum. 
Since spectral efficiency of optimum joint decoding is an increasing function of the load [4], anti-resource 
pooling is beneficial for spectral efficiency, though its implementation may cause some practical challenges. 

In the following theorem, we extend anti-resource pooling to arbitrary delay distributions: 

Theorem 2 Let A G C fj e a diagonal matrix with k th diagonal element G C and T c a positive 
real. Given a function $(a>) : R — >• C, let 4>{VL,t) be as in 4721) . Given a positive integer r, let <fr k , 
k = 1, . . . , K, be r -block-wise circulant matrices of order N defined in ([77]). Let H = SA with S = 
[<&isi, $2^2, . . . , <&kS K ] and s k G C Nxl . 

Assume that the function \<&(u)\ is upper bounded and has support contained in the interval 
The receive filter is such that the sampled discrete-time noise process is white. The vectors s k are indepen- 
dent with i.i.d. circularly symmetric Gaussian elements. Furthermore, the elements a k of the matrix A are 
uniformly bounded for any K . The sequence of the empirical distributions F^l(X) = ^Yl^=iX(^ — l a fc| 2 ) 
converges in law almost surely, as K — > oo, to a non-random distribution function -Fj^^A). 

Then, the multiuser efficiency of the linear MMSE detector for CDMA with transfer matrix H converges 
in probability as K, N — > oo with £ — * j3 and r fixed to 

+tt/T c 

lim f] k = rj = — / 7] (a>) dcu (35) 



7T 7T 



K=PN^oo 2lT 

-tt/T c 

where the multiuser efficiency spectral density rj(uj) is the unique solution to the fixed point equation 

1 _ E, + £f^WW (36) 



n{u) |*M| 2 tJ % + 

for all uj in the support of^{u) and zero elsewhere 
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Theorem [2] is proven in Appendix ITVl 

No constraint is imposed on the set of time delays in Theorem[2l It holds for any set {t±, r 2 . . . r K } and 
we conclude that linear MMSE detectors for synchronous and asynchronous CDMA systems have the same 
performance if the bandwidth of the chip pulse waveforms satisfies the constraint B < ^r. 

VI. Spectral Efficiency 

There exists a close relation between the total capacity of a CDMA system and the multiuser efficiency of a 
linear MMSE detector for the same system 0, 11361 . [|37l . The rationale behind this relation is a fundamental 
connection between mutual information and minimum mean-squared error in Gaussian channels ll39ll . In the 
following, we extend the results in Section |V] to get insight into the spectral efficiency of an asynchronous 
CDMA system. 

The capacity of the CDMA channel was found in [40] for synchronous CDMA systems. The total capacity 
per chip for large synchronous CDMA systems with square root Nyquist pulses and random spreading in the 
presence of AWGN (additive white Gaussian noise) is H| 



C (syn) (/3,SNR) =/51og 2 (\ + SNR-^F(SNR,/^ + log 2 ^1 + /?SNR - ^F(SNR, /3)J 

-3i§g«> (37) 

with 

F(y,z)= (y^l + v^ + l- yVl-v^ + l) • (38) 

With the normalizations adopted in the system model, we have SNR = E^/Nq. 

The spectral efficiency of a synchronous CDMA system is equal to C^ syn ^(/3, SNR) for any Nyquist sine 
waveform. For other chip waveforms, we need to take into account the excess bandwidth and calculate 
spectral efficiency as 

r = h m 

where C denotes the total capacity per chip and B denotes the bandwidth of the chip pulse. Note that for 
Nyquist sine pulses T C B = 1, while in general T C B can be either larger, e.g., for root-raised cosine pulses, 
or smaller, i.e., for anti-resource pooling, than 1. 

The expression of the total capacity per chip for asynchronous CDMA systems constrained to a given chip 
pulse waveform ip(t) of bandwidth B and a given receive filter g(t) can be obtained by making use of the 
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results in Section |V] and the fundamental relation between mutual information and MMSE in Gaussian chan- 
nels provided in ll39ll . Since such constrained total capacity depends on ip{t) and g(t) only via the waveform 
output of the filter g{t) for the input ip(t), we shortly refer to it as the total capacity constrained to the 
chip waveform 

Corollary 4 Let us adopt the same definitions as in Theorem [2] and let the assumptions of Corollary [2] 
or Theorem [2] be satisfied. Additionally, let the receive filter and sampling process be such that sufficient 
discrete-time statistics are provided. Then, as K, N — > oo with ^ — > {3 the total capacity per chip con- 
strained to the chip pulse waveform (fi(t) converges to the deterministic value 

o 

where r/ 7 is the multiuser efficiency at signal-to-noise ratio 7 given in rf25l) and 051) . respectively. 

The proof of this corollary is discussed in Appendix IVl! 

Let us consider again the case of sine chip waveforms as defined in (|3TT) and uniform distribution of the 
time delays. Let a denote the bandwidth of the sine pulse relative to the Nyquist bandwidth. As noticed in 
Section |V] the multiuser efficiency ?7 S i nc of an asynchronous system with such sine waveforms given by (1331) 
and load (3 equals the multiuser efficiency ^ syn of a synchronous system with Nyquist sine pulses given by 
(|34l) and load /?' = -. Since the load enters capacity per chip (|40b only via the multiuser efficiency except for 
the linear pre-factor to the integral, we immediately find the following equation relating the two capacities 
per chip 

C (sinc) (/3, SNR, a) = «C (syn) (?-, SNr) . (41) 

It is apparent from (PTTI) that synchronous and asynchronous systems have the same capacity for a = 1. 

In order to compare different systems (with possibly different spreading gains and data rates), spectral 
efficiency has to be given as a function of the level of energy per bit per noise level equal to H| 

E> = (3SNR 

N CM (A SNR,-)' 

In Figure [2l we compare the spectral efficiency of asynchronous CDMA with the spectral efficiency of 
synchronous CDMA. TheF spectral efficiencies are plotted against the bandwidth normalized to the Nyquist 
bandwidth with = 10 dB and unit load (3 = 1. Recall from earlier discussions that for synchronous 
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Fig. 3 



Spectral efficiency of random CDMA with unit load Spectral efficiency of random CDMA versus the load p 
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No 

STANDARD AND = 10 DB. 



systems all Nyquist chip waveforms perform identically. So there is no need to specify a particular Nyquist 
pulse except for the Nyquist pulse having the same bandwidth than the sine pulse in the asynchronous case. 
We see further that the smaller the normalized bandwidth, the higher the spectral efficiency is. This is, as 
anti-resource pooling improves spectral efficiency by emulating a higher load. 

In Figure [3] the spectral efficiency is plotted against the load (3 with jj^ = 10 dB for the chip waveform 
used in the UMTS standard. When the load (3 increases the gap in spectral efficiency between synchronous 
and asynchronous systems increases. 



VII. Extension to General Asynchronous CDMA Systems 

In this section we extend the previous results to any distribution of the time delays for CDMA systems. 
Without loss of generality we can assume that the time delays r k E [0, T s ) [|28l . In this case, intersymbol 
interference is not negligible and an infinite observation window is necessary to obtain sufficient statistics. 
Equation (fl4l) for the chip asynchronous but symbol quasi-synchronous system model is extended to a gen- 
eral asynchronous system by 
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K +00 



Tk 



w[p] 



(43) 



with p E Z and 



fc=i 



4 m) = E s *,™M# - (^+^)t c ) 



JV-1 



(44) 



u=0 



By assuming the same approximation as in (fl4l) . the virtual spreading sequence of user k in the symbol 
interval m has nonzero elements only in the time interval m and m+1. Let Tk denote the delay of the signal k 



and r k = r fc modT c , respectively. 



in terms of the chip intervals and r k the delay within a chip, i.e., Tk = 
The virtual spreading sequence of user k is obtained by computing <3? fc as in (fTTI) for r = r fc to account for 
the delay within a chip and then by shifting the virtual spreading vector down by Tk r-dimensional blocks to 
account for the delay multiple of the chip interval. More precisely, the virtual spreading in the m-th symbol 
interval is given by the 2r iV-dimensional vector 



with 



C[T k ) 



0T k 



,M 



(45) 



T, 



(r - 1)T C 



0- k and 07v_- fe column vectors with zero entries and dimension rr k and r(N — r k ), respectively. The 
2r N x K virtual spreading matrix for the symbols transmitted at time interval m is then 



;(m) 



— (m) ~ (m) ~ (m) 

For further study, we introduce the upper and lower part of the matrix S , S u and S d of size r iV x rK 



such that 



dm) 



dm) 



dm) 



{in) ~(m) (m) ~(m) 

and the matrices i/ u = S u A and ff d = S d A. Then, the baseband discrete-time asynchronous 
system in matrix notation is given by 



y = TiB + W 



(46) 
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where y 



[■■■,y 



(m-l) T 7| (m) T ,.(m+l) 



y 



...fandB = [. . . M m ~ lY 'M m)J \b 



(m.+l) 1 



} T are the infinite- 



length vectors of received and transmitted symbols respectively; W is an infinite-length white Gaussian 
noise vector; and 7i is a bi-diagonal block matrix with infinite block rows and block columns 



n 



o // 



(m-l) 



H 



(m) 



H 



(m) 



H 





(m+l) 







(47) 



Finally, we define the correlation matrices T = TVH H , 71 = H H H. 

The following theorem shows that a linear MMSE detector for a CDMA system with transfer matrix 
7i and time delays r±, T2, . . .tk nas the same limiting performance as a linear MMSE detector for chip 
asynchronous but symbol quasi-synchronous CDMA systems introduced in Section [IV] with time delays 
ri, r 2 , . . . tk- The same equivalence holds for capacity and spectral efficiency. 



Theorem 3 Given {n, r 2 , . . . tk} a set of delays in [0, T s ) let us consider the set of delays in [0, T c ) de- 
fined as {% : % = Tk mod T c , k = 1, . . . K} . Given a positive integer r, let k = 1, . . . K, be the 
r-blockwise circulant matrix of order N defined in 4771) with r = r fc . Let A, $(cj), S, and H be defined 
as in Theorem\l\ Furthermore, k — 1 ... K are 2rN x N matrices such that = [0^ , 0]y_-J T 



with Tk = 
Let S im) 



Ik. 



0r k and 0N-r k zero matrices of dimensions rTk x N and r (N — Tk) x JV, respectively. 



(m) 



,(m) 



~(m) r ~(m)T ~(m)T. 



S A an J ?i infinite block row 



and block column matrix of the same form as in fflTb . Let the same assumptions as in Theorem{T}hold. 

Then, asymptotically, as K, N — > 00 with jt — > j3 the CDMA systems transfer matrices Ti and H are 
equivalent in terms of multiuser efficiency for linear MMSE detectors and in terms of spectral efficiency. 



This theorem is shown in Appendix I VIII 

Interestingly, the system performance depends on the time delays r fe only through the offsets r fc — 
Therefore, any shift of the signal multiple of T c does not affect the performance of the system. 

The analysis presented in this contribution has been restricted to frequency flat fading for the sake of 
clarity. The extension to multipath fading channels is straightforward when the impulse response of the 
channel is much shorter than the symbol interval. In fact, the chip pulse waveform at the output of the 
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matched filter 0(t) can include the effects of the frequency selective channel impulse response a(t) along 
with the effects of the transmitted chip pulse waveform ip(t), and the filter at the front-end g(t), i.e., 4>{t) = 
ip(t) * a(t) * g(t). Then, the analysis of a system with frequency selective fading reduces to the proposed 
analysis. 

VIII. Conclusions 

This work provides a general framework for the analysis of asynchronous CDMA systems with random 
spreading using sufficient or suboptimum statistics and any chip pulse waveform. Furthermore, it includes 
several optimum or suboptimum receiver structures of practical and theoretical interest. Therefore, it pro- 
vides insight into both the fundamental limits of asynchronous CDMA systems and the performance loss of 
implementations where suboptimum receiver structures, suboptimum statistics, and/or non-ideal chip pulses 
are utilized. 

For the receiver structures investigated in Part I, the performance of a CDMA system is independent of 
the time delay distribution if the bandwidth of the chip pulse waveform is not greater than half of the chip 
rate, i.e., B < This also implies that synchronous and asynchronous CDMA systems have the same 
performance and generalizes the equivalence result in lfT8l for Nyquist sine (B = ^r) pulses and linear 
MMSE detectors to any chip pulse waveform. The behavior of CDMA system changes substantially as 
the bandwidth gets larger. In this case, the system performance is significantly affected by the distribution 
of the time delays and the performance of linear detectors may depend on the specific time delay of the 
signal of interest. If the receiver is fed by sufficient statistics and the time delay distribution is uniform the 
performance of optimum or suboptimum receivers is independent of the time delays. In the following, we 
summarize the most interesting aspects pointed out by the large system analysis, for each class of receivers. 

A. Optimum Receiver 

The spectral efficiency constrained to a given chip pulse waveform characterizes the performance of a 
CDMA channel with optimum receiver. The spectral efficiency is expressed in terms of the multiuser effi- 
ciency spectral density t]{oj). When the chip-modulation is based on sine pulses whose bandwidth is a times 
the Nyquist bandwidth, the spectral efficiency of asynchronous CDMA systems is identical to the spectral 
efficiency of synchronous systems with load f3' = - and Nyquist sine pulses. Spectral efficiency is a strictly 
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decreasing function of the relative pulse bandwidth a and for a — > 0, the spectral efficiency of a single user 
AWGN channel is reached. 

For a > 1 an asynchronous CDMA system with modulation based on a sine function can compensate to 
some extent for the loss in spectral efficiency of synchronous CDMA systems with equal bandwidth. For 
(3 — > oo it attains the maximum spectral efficiency for any finite bandwidth B = ^fr. 

B. Linear MMSE Detector 

The output SINR of a linear MMSE detector can be obtained from the solution to a system of fixed point 
equations in the general case. In the two cases (i) chip pulses with bandwidth B < and (ii) chip pulses 
with bandwidth B > sufficient statistics and uniform time delay distribution the fixed point system 
of equation reduces to a single equation. In those cases, the performance of a linear MMSE detector in 
asynchronous CDMA systems is characterized by a unique value of multiuser efficiency. Furthermore, the 
measure of multiuser efficiency can be refined by the concept of spectrum of the multiuser efficiency that is 
also unique for all the users. Additionally, for these CDMA systems the limiting interference effects can be 
decoupled into user domain and frequency domain, as the system grows large, such that we can define an 
effective interference spectral density similarly to the effective interference in @ for synchronous systems. 

In the special case that the modulation is based on sine functions with bandwidth B — a linear 
MMSE detector in asynchronous CDMA channels performs identically to a synchronous CDMA system 
with square root Nyquist chip pulses J51 and load f3' = -. This effect is similar to the resource pooling effect 
for synchronous CDMA systems with spatial diversity in |[38Tl and shows the possibility to trade degrees of 
freedom in the frequency domain against degrees of freedom in the time domain. 

Though this work focused on performance measures for CDMA, similar results hold for asynchronous 
MIMO systems due to the mathematical analogy between CDMA and MEMO systems when described as a 
discrete-time vector channel. This means, that MIMO systems with excess bandwidth and desynchronized 
modulators for different antenna elements benefit in a similar manner than CDMA systems with desynchro- 
nized users. 
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Appendix I 
Useful Mathematical Tools 

Let $(o;) be the unitary Fourier transform of a pulse waveform <p{t) with bandwidth B < ^r. Then, in 
the normalized frequency interval Vt E [— ir, n] the unitary Fourier transform (fL2l) of the sequence obtained 
by sampling 4>(t) at time instant r and rate ^- is given by 



*(!),r)=ie«« J) ^**(^) 



for Q < 7r 



(48) 



where Z{VL) is the set of all integers in the interval 

[min (-sign(fi) , sign(fi) [fj) , max (-sign(fi) [^J ,agn(fi) [|J)] . 

The matrix 

Q(n,r) = A^.(n,r)A^ r (n,T) H , 
with A^ r (fi, r) defined in (1221) . can be decomposed in the sum of two matrices 

Q(n,r) = Q{Q) + Q(Q,t) 

where the (k, £)-elements of the matrices Q(fi) and Q(£l, r) are given by 



(49) 



(50) 



se2(n) 



f2 + 27TS 



Tr 



-j^(n+2ns) f()r , n . < 



7T, 



(51) 



and 



n + 2nu\ _ / r fi + 27ra\ e _ J .^ (a _ B)e _ i (fc=i (n _ 2w , ) _fci (n _ 27ro) ) 



s,uS.2(Q) 



for Ifil < 7T, (52) 



respectively. 

Useful properties of the matrices Q(fi) and t) are stated in the following lemmas. 



Lemma 1 Let B be an r x r matrix of the form 



B = B(Q) 



-ie 



b b\e^^ 



(.r-l) , 



6 r _ie" 



(53) 



November 28, 2009 



Submitted manuscript to IEEE Transactions on Information Theory 27 
i.e., given b = bo(Q), b\ = b\ (fi), . . . 6 r _i = 6 r _i(f2), eventually functions ofVL, (B)i^, the element (£, k) 
of the matrix B satisfies (B)^ = e ? n b( r+ k-e)modr- Let Q(Q, r) be the r x r matrix with element (k, £) 



defined in 021) . Then, 

ti(BQ(n,r)) = 0. 
Proof: Letq us (n) = ^$ $* (^r) ' Then ' 

r r 

tv(BQ(n + j27ru, r)) = E^, T)) k/ (B) e , k 

k=l 1=1 



7^(-«(*-l)+«(<-l)) 



s,ue.z(n) fc=i f=i 

E ?« s e^(*^)^«) E E 6(r+M)-»dre-^ ( "-V^- 



with 



and 



E ?«y 27r( *^ )(s_w) (^+^) 

s,«G.2(f2) 



Vl = b (r+k-£)mod r e rK ' Z r K ' (54) 

k,e=i 

k>l 



V 1 h a -3 — {u-s)kJ^-u(k-e) 



k,e=i 

k<t 

Substituting v = k — £ in (l54l) and v = r + k — £'m (l55l) we obtain 

r fc-1 



-j22L(u-s)fe Q j2zni£ 



k=l u=0 

and 



fe=i u=fe 



respectively. For s,i 6 and s 7^ i, |s — t| G [1, . . . , r — 1]. Therefore, Y7k=i e~ J V"( u ~ s ) fe = and 

^1 + = for all f2. Then, also tr(I3Q(fi, r)) = and this concludes the proof of Lemma[TJ ■ 
It follows immediately from Lemma \T\ that trQ(Q, r) = since the identity matrix I is of the form (|53l) 
with b — 1 and 6j = for % — 1, . . . r — 1. 
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Lemma 2 Let B = B(Q) be a matrix defined as in LemmaU\and let Q(Q) be the r x r matrix with element 
(k, £) defined in ( 1571) . Then, the matrix C(Q) = Q(Q)B(Q) is of the form rt53l) . 

Proof: The element (k, £) of the matrix C = C(fi), (C)^ is given by 

r 



)t,k 



t=l 



e r 



k(£, k) 



(56) 



with 



and 



sez(n) 



r}(£, k,s) 



T}(£, k,s)=^ b (r 



(57) 



t=i 



In order to prove Lemma |2] it is sufficient to prove that 



k(£, k) = k((£ + l)modr, {k + l)modr). 



(58) 



In fact, in this case C^ k = e> Q K( r +k-e)modr with K {r+k-e)modr = k(£, k). The property (|58l) is implied by 
a similar property on r]{t, k, s) 



r)(£, k, s) = r)((£ + l)modr, (k + l)modr, s). 



(59) 



It is straightforward to verify that ([591 is satisfied since both factors 6( r +t-^) mo dr an d e^M ~) are periodical 
in their arguments £ and k, respectively, with period r and k and £ are simultaneously increased by a unit. 
This concludes the proof of Lemma[2] ■ 

The following lemma provides the eigenvalue decomposition of the matrix Q(fi). 

Lemma 3 Let Q(Q) be an r x r matrix with element (k, £) defined in ( 1571) . Then, the matrix Q(Q) can be 
decomposed as follows 

Q(n) = u(n)D(n)u H (n) (60) 



where 



U(Q) =(e(n- sign(fi)2vr 



r — 1 



,...e(ft)...e m + sign(f2)27r 



(61) 
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e (Q) is an r- dimensional column vector defined by 



and D(Q) is the diagonal matrix whose s th diagonal element is given by 



[Dm 



r 

SS rj-,2 



$ I y ( Q - sign(fi)27r 



r — 1 



-s + 1 



(62) 



Proof: Decomposition (1601) can be immediately derived by noting that 

2 



r 

J^2 



I 7^ + 27TS 



Q(0)= £ 

This expression can be rewritten as (l60l) and Lemma[3]is proven 



e(n + 2irs)e H (n + 2irs). 



The following lemma shows that the matrix Q(fi) and any other matrix with the same basis of eigenvectors 
is of the form (|53l) . 

Lemma 4 Le? C(fi) = U(Q)M(Q)U H (Q) with U(Q) unitary matrix defined in rfoTT ) and M(fi) diagonal 
matrix with elements m/ c / (f2). Then, C(Q) is of the form ( T53l) . 
Proof: The £ th row of the matrix £/(f2) is given by 

= _L ^ e -i £ T i (n-sign(n)2^L 2: i i j) ) _ _ _ e -j^(n+sign(n)2^L|j) 

and Qfc(fi), the element (£, k) of the matrix C satisfies 

5 (fi) = I^ m . i e-^( n - si ^) 2 -L^J+ 2 ^- 1 )) 



Qfe 



i=l 



~ 7 -(*-»0o 



(63) 

)modr • 



with b lk = EI=i ^ e J ' 2 ^( sign ^^J-+i)_ It is straightforward to verify that b ek = b {e+1)modr , ik+1 
This concludes the proof of Lemma|4] ■ 
The following lemmas state results from random matrix theory developed along the lines of the REFORM 
method proposed by Girko in iHTI and B2l . 

Lemma 5 K4T] . K43\l Let E = (Cii)i=i'.'.'.JVgi ^ e an x matrix of complex random elements &j 
structured in NK blocks of size qi x q 2 , 3y, z'.e., 

w _ cw..y=i.-^ 
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andK = (3Nwithf3 > O.LetP = (Ptj)y=i,...pi = [HH^+aJ]^ 1 andG = (C? i j)ii=i ) ...p a = [S^H+a/]" 1 , 
where P^ and Gij are complex blocks of size q\ x q 1 and g 2 x Q2, respectively. 
Additionally, assume 

H-l Sfes, fc = 1, . . . , N, s = 1, . . . , K, the random blocks of the matrix E are independent. 
H-2 All the elements of the matrix E are zero mean, i.e., E{E} = 0. 
H-3 sup^max i= i ) ... )J vI]jL 1 E||S ii || 2 + sup^max^i,...^ ^JIi E ll 3 iill 2 < +°°> 
H-4 Lindeberg condition: Vr > 

kJ^oo y i ^ N Jl E (\\^\\ 2 ^\\ S i3\\ > r }) + j ^ K J2 E (HSijIlMllSyll > r}) J = 0. 

(64) 

Then, for a G C\R _ 

lim E|P p£ (a)-T p£ («)| = p^=l,..., Pl 

K=pN — >oo 

lim E|G P £(«) - a~ l R pl (a)\ = p,£ = l,...,p 2 

K=f3N^oo 

i.e., the blocks of the matrices Q and G converge in the first mean to the corresponding blocks of the matrices 

f = diag((C«(a))- 1 ) n=1 ,... iAr 

and 

R = dmg((C%(a)r 1 ) k=h ...,K 
respectively. The matrix blocks C£2(a) of size q\ x q 1 and C^(a) of size g 2 x q2 are equal to 

K 

C<2(a) = oT + E ( S ^( X )«' S S) X=Q 6 n = 1, . . . , JV (65) 

j'=i 
pi 

CS (a) = I + X) E (Hf fc (^),,Hf fe ) y= - fc = 1, . . . , K, (66) 

3=1 
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Lemma 6 / l?71/ . I[43\l Let us assume that the definitions of Lemma\5\hold and the conditions of Lemma\5\are 
satisfied. 

Then, the qy x q x matrices C^(a), n — 1, . . . , N and the q 2 x q 2 matrices C^(a), k — 1, . . . , K, defined 
in ( 1371) and ft66\) , respectively, converge as K = j3N — > oo to the limit matrices 

lim CS=#Cg n = l,...,JV 
lim c2=*2 fc = l,...,A: 

K=/3N->+oc KK KK 

where VE*^, A; = 1, . . . , N and VE^ , = 1, . . . , if satisfy the canonical system of equations 



V^al+fy^nj [ttg^S*} , n = l 



A' 



iV, (67) 



-i 



*i?=^+E E S SKJ ' = (68) 

The following Lemma states the existence and uniqueness of the solution of the system of canonical equa- 
tions in the class of definite positive Hermitian matrices. 

Lemma 7 K41\l Let us adopt the definitions of Lemma |5j and let us assume that the conditions of Lemma 
UJ are satisfied. Let us consider the system of canonical equations WJ\ and (US]). Then, the solution of 
the canonical system of equations ft67\) and ft6E\) exists and it is unique in the class of nonnegative definite 
analytic matrices for Re(a) > 0. 

The following lemma due to Girko provides convergence of the eigenvalue distribution of the matrix HS^ 
with S defined in Lemma[5]to a deterministic distribution function and the corresponding Stieltjes transform. 



Lemma 8 /HZ]/ Let us adopt the definitions in Lemma\5\and let the assumptions of Lemma\5\hold. Further- 
more, let n qi N(x, HS^) denote the normalized spectral function of the square q\N x q\N matrix argument, 
i.e., the empirical eigenvalue distribution of the matrix HS^. Then, for almost all x with probability one, 

lim \n qiN (x, EH^) - F qiN (x)\ = 

N— >oo 
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where F Qi n(x) is the distribution function whose Stieltjes transform is equal to 

r+oo ^ 

/ (x + ay'dF^ix) = (q 1 N)-Hv[^\- 1 (69) 
Jo 

with = diag(^ r nn ) n=1 jv nonnegative definite analytic matrix for Re(a) > and ^ nn satisfying the 
canonical system of equations rfoTl ) and rfolSl ). 

Lemma 9 Let x = (xi,X2, ■ ■ ■ , x^) be an N -dimensional column vector of complex i.i.d. elements 
with zero mean and unit variance and C be an N x N complex matrix. Then, for any p > 2 

E\x H Cx - trC\ p < K p ((E|x!| 4 trCC^) f + (E|:n| 2p tr(CC H ) (70) 

with K p positive constant independent ofN. 

Appendix II 
Proof of Theorem CD 

Let us consider the r-block-wise circulant matrices of order N, k = 1, . . . K defined in Theorem [H 
and let us denote with F 1 ^ the unitary Fourier transform matrix of dimensions N x N with (I, m) element 
given by 



(F N ) e , m = ^=e^^l (71) 
v N 



rices! [1311 



We can extend the well known results on the diagonalization of circulant matricesa DTI to decompose the 
r-block-wise circulant matrices <3> fc , k = 1, ... K as 

= {F N <g> I r )A^ r {r k )F N (72) 

where A^^rfe) is an rN x N block diagonal matrix with £ th block given by 

(A , r (r fc ))^ = (svr^i, r fc ^ (73) 

and (Fjy ® I r ) is a unitary matrix. 
The matrix S can then be rewritten as 

S = (F N <g> I r )(A^ r (ri)5i, A0 )r (r 2 )s 2 , • • • , A^ r (T K )s K ), 

5 A circulant matrix C = 3(f(x)) of order iV can be decomposed as C = F N DF%, with D = diag(/(0), /(^), . . . , f(2n^^-)). 
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with s k = F^Sk. Assuming the elements of the spreading sequence s k i.i.d. Gaussian distributed, s& is also 
a vector with i.i.d. Gaussian distributed elements having the same distribution as the elements of s k . Since 
the eigenvalues of any matrix X are invariant with respect to left multiplication by a unitary matrix U and 
right multiplication by U H , i.e., the eigenvalues of the matrix X coincides with the eigenvalues of the matrix 



UXU H , then the singular values of the matrices S and S = (A^, r (ri)si, A0 jr (r 2 )s2, 



'Jk)s k ) 



coincide. The same properties holds for the matrices H and H = SA. It is straightforward to verify that 
also SINRfc is invariant with respect to such a transform. In fact, 

SINR fc = hf (H k H? + a 2 !)' 1 h k 



a k \ 2 s^[H k H k +0*1 



H 



-1 



with H k and H k obtained from the matrices H and H , respectively, by suppressing the k th column. There- 
fore, in the following we focus on the analysis of the system with transfer matrix H. 

The matrix H is a matrix structured in blocks of dimensions r x 1. The block (n, k) h n k, n = 1, . . . N 
and k = 1 , . . . K , is given by 

h"n,k \^ j k\ (.^(j}, r (.^~k)^)rinSn,k 

where J Thk is a Gaussian random variable with zero mean and variance Eds^] 2 } = 4. Additionally, the 
variables s n fc are i.i.d.. Therefore, conditions H-l and H-2 for the applicability of Lemma[5]and Lemma[6] 
are satisfied. Condition H-3 of Lemma[5]is satisfied. In fact, 

K N 

max VE{||h„ fc || 2 }+ max VE{||h nfc || 2 } 



C = sup 

N 



< sup 

JV 



fc=l 

K 



71=1 



K | , 2 | , 2 JV 

m ^ T ^-^\\i A 4>,r{h))nn\\ 2 + max V || (A^ r (r k )) nn 

n=l...N ^— ' iv k=l...K iv 



fe=l 



n=l 



Since the function $(c<j) is bounded in absolute value with finite support also |$(fi, r) | is upper bounded for 
any f2 and r. Then, there exists a constant Cmax > that satisfies || (A0 jr (r fe )) n „|| 2 < Cmax for any k and 
n. Additionally, the elements \a k \ are uniformly bounded for any k, i.e., 3omax > such that \a k \ 2 < a MAX 
for all k. Then, 

C < sup a MAX C M AX ( ^ + 1 J < +oo. (74) 

X=/3JV V ^ J 

In order to verify the Lindeberg condition H-4 we focus on the limit 



K=(3N 



K 

l ^^ m ^J2 E (^k\\\(\\Kk\\ >*)) 
fc=l 
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for any 5 > 0. Let us observe that Vn, k 

E ( ||^nfc|| 2 X(||^nfe|| >£)) = |afc| 2 ||(A0 jr .(r fc )) nrl || 2 / ISn/^dF^fc) 



k| 4 ||(A^(r fc )) I' 4 



< / \s n k\ dF(s nk ) 



d 2 



{\s nk \>0} 



where F(s nk ) is the cumulative distribution function of s nk . By using the fact that s nk is a complex Gaussian 
variable with variance Eds^ 2 } = and forth moment E{|s n fc| 4 } = ^ and the bounds on |a fc | 2 and 



'A^ r (7fc)) nn ||,it holds 



[\\h nk \\\{\\h nk \\ > 5)) < 2a MAAx _ (y5) 



Then, 7] = since 



A' 



o < 77 < 2q4ma ^ax Hm ri 

52 A-=/3JV->oc ^ TV 2 

fe=l 



Similarly, it can be shown that 

lim max V E ( \\h nk \\ 2 x{\\h n k\\ 2 > 5)) = 

K=f3N^oo k=l,...,K ^ V / 
n=l 

and the Lindeberg condition H-4 is satisfied. 

From Lemma[51 U p e(a), p,£= 1, . . . , N, the blocks of the matrix U(a) = (H k H k + al)^ 1 converge in 
the first mean to r x r matrices V p g = (Cf e ) )~ 1 5 p e, p, £ = 1, . . . N, and C$ defined similarly as in Lemma 
[5] Additionally, from Lemma [6] the matrices C^} can be obtained as solution of the canonical system of 
equations (1671) and (1681) asymptotically as K = j3N — > 00. Equations (1671) can be rewritten as 

K 

Y« = a J r + ^E{^ fc [Ti 2) ]- 1 v} 

fc=i 

= a/ - + ^ D T 2] _1 l«*l 2 **r ( 27r! ^r ' r *) A £ ( 27r! ^r ' " = L • • • N ™ 

with (x, r) defined in (1221) and taking into account that (A^ r (r fc ))„„ = A^ r (27r 2 ^-,r fe ) in (|76|) . 
Equations (1681) specialize to 



T (2) — 1 + 

n=l 

I I9 JV 

= 1 + 



g Aj r (2^^, r fc ) [T«]"' A,, r (2^^, t*) k = 1, . . . K. (77) 
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By substituting (1771) in (l76l) and considering the canonical system of equations as K = (3N — > oo we obtain 

A (fi, r) A l r (Q, t) f\ A? , T {\ r)d Ad r 



m /" AAj, - [it, T) A 



w J* A J r (tf, r) [T (1) (fy)]"* A Ar (tf, r) d fi' 
with G [0, 2ir], or since A<^ iT .(f2, r) is periodical in f2 with period 2n, f2 can equivalently varies in the 
interval X = Here, S denotes the support of the distribution function F\ A \2 T (X, r). By defining 

T(fi) = [T (1) (fi)] _1 we obtain (ED). It follows from Lemma[6] 



K=PN 



lim Cffi = T 1 (2^) . 



-H 



The convergence in the first mean and thus in probability of SINR^ = h k U(o )h k to the quantity g 

-H 

h 



^ / A£ r (fi,7*)T(fi)A , r (fi,7*)dfi is proven if Vl = E 



KU{a 2 )h k - g 



vanishes asymptotically, i.e., 



lim 77i = 0. 

jv ^ 



(78) 



The rest of the proof is focused on showing (1781) . Let us observe 

77! <E 



h H h U(a 2 )h k - h!.Vh k 



~H ~ 

E|h fe Vh k — g\ 



(79) 



"k ^ \ u - '"k 

where the triangular inequality^ of the spectral norm is applied and V = diagQC&V 2 )]- 1 )* =l,...,jV 1S 
defined in Lemma [6] 

By applying the submultiplicative inequality for spectral norms and the triangular inequality to the first 
term of (|79l we obtain 



'H 



Eh t I/ O - V)h k \ = E 



i.t 



^kKr ( 2vr^,r fe ) (U(a 2 ) - V) u A^ r ( 27r^,r fc ) 



N 



N 



< J2 E (\(U{v 2 ))it- Vu\)E\r ik se k \A^ r (2n l -^,r k ) A^ r (2/-^,, 



N 



\A^r (2lT^,T k )\ 

N 



< E 11 ^" (27r /' Tfc) l|2 maxE|(£/(a 2 )), - V u \. 
z — ' iV t 

6 Given two matrices A and B with consistent dimensions the following inequalities hold: 

\AB\ < Submultiplicative inequality of spectral norms; 

A + B| < |A| + |2?| Triangular inequality of spectral norms. 



N 
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Thanks to Lemma[5]and the fact that || A^ (2ii^-, r fc ) || 2 < K M ax for all z = 1, . . . iV and r k 

lim E\h H k (U(a 2 )-V)h k \ = Q. 

K,N— >oo 

n y 

~H ~ 

In order to prove the convergence to zero of rj 2 = E|/i fe Vh k — g\ we consider 

V 2 < E|hf Vh k - g\ 2 

= E((h"vh k ) 2 - 2gh"vh k + q 2 ) 



= E ( M 4 £ A*. (27r^,T fe j V«A, )P (2vr^,r fc J A* r ^tt^^J V„A.,, ( 27 r^,r fc ) |^| 2 M 
-2g|a fe | 2 ^A^ ^27r^,^ V«A , r ^tt^t^ |^ fe | 2 + g 2 j (80) 

X ^ ^V' Tfc ) (^^'^ - ( 2 ^'^) ( 27r T' Tfe ) + ^ (81) 

From (|80l ) to (18TI ) we make use of the assumption that 3^ is a complex Gaussian variable circularly invariant 
with variance N^ 1 . Let us observe that the spectral norm of T(f2) and Va, for any i, are bounded by 
|T(f2)| < a 2 and \ Va\ < a 2 . Then, the first term in (I8TT) vanishes as N — > oo. By appealing Lemma[6l for 
any i, — > T (27r^) as .fT, JV — >• oo with — > (3. Then, the second and third terms in (|81"I) converge to 
£ 2 and — 2g 2 1 respectively. We can conclude that 

lim ?i 2 = 

K,N^oo 
N " 

and ?7 2 — > as if, N — > oo as — > /3. Therefore, (|78l) and thus the convergence in the first mean of SINR fc 
is proven. The Markov inequality implies that, Ve > 

lim Pr {\h*U (a 2 ) h k - g\ > e} < - lim E|hf U(a 2 )h k - g\ = 

K,N—>oo £ K,N—>oo 

^^0 

and the convergence in probability stated in Theorem Q] is proven. 
This concludes the proof of Theorem [Q 
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lim SINR fc = ^/ tr(T(ft)Q(fi,7fc))dft (82) 
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Appendix III 
Proof of Corollary |2] 
In order to prove Corollary [2] we rewrite the limit SINR^ in (1201) as 

3 INR t = 

and the fixed point equation (l2"TT) as 

Jo Ja l + ^/!,tr(r(ti')Q(!i',T))dti' - - 

with Q(fi, r) defined in ( |49l ). The matrix r) can be decomposed as in ( |50l ). Thanks to the assumptions 
on in Corollary [2l the conditions on Q(f2) and Q(f2, r) in Lemma [Hand Lemma [2] are satisfied. 

First we show that Y(fi), the unique solution of (1831) in the class of nonnegative definite analytic functions 
in Re(<7 2 ) > 0, is an r x r matrix with eigenbasis U(fi) defined in (loTI) . Let us assume that Y(fi) = 
C/(fi)T(fi)?7^(fi) with elements of T(fi) nonnegative for all Q e [— 7r,7r] . By appealing to Lemma |4] 
Y(f2) is of form (153T ). Then, by applying Lemma[l]it results tr (Y(fJ)Q(fi, r)) = for all f2 G [— 7T, it] . 
Therefore, 



(83) 



/7T 
tr(r(tt')Q(tt'))dtt' 
■n 

tr (Y(fi / )£>(^')) dfi/ > 



with D(Q) defined as in Lemma[3l Let us notice that J^ c Q(Q, r)dFy(r) = for all fl Thanks to this 
property, the assumption of independence of the random variables A and r, and to the uniform distribution 
of r, (1831) can be rewritten as 

Y (SI) = a 2 I r +P\ ' 1 1 ; r - 7T < ft < TT. (84) 




Since Y(ft) and are diagonal matrices, the matrix equation (1831) reduces to a system of L scalar equa- 
tions. Furthermore, all the quantities that appears in the right hand side of the system of equations < [8"4l ) are 
nonnegative under the assumption that Y(ft) is a nonnegative definite matrix and (l84l) admits a nonnegative 
definite solution for Re(<7 2 ) > 0. The existence of a nonnegative definite solution of the system of equations 
(l84l) implies also a solution of the fixed point matrix equation (f83l) given by Y(ft) = [/(ft) Y (ft) t/^ (ft). 
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Let v s (n) be the s diagonal element of T(fi) and let us recall that the s diagonal element of D(Q) is 
given in ( 1621) . Then, < [84l ) reduces to 



$ ( — -sign(fi) — 





r — 1 


-0) 


( 






2 





X 



-7T < < vrands = l,...r. (85) 

By changing the variable y = f2 — sign(f2)27r ( — s + l) and defining the function v(y) in the interval 
(— m, m) as follows 



v(y) 



v. 



(y-2n([^\-s + l)) ~s + \<l< L^J - « + 



".(v + MI/^J "* + !)) S -l-Lr=lj<JL< a _l_Lr=lj 



the system of equations (|85l ) can be rewritten as 

2 -+00 



" 2 + % 



A similar approach applied to (1821) yields 



Ad J F| Ap (A) 



Urn SINR fc = ^^ 

K,N->oo T? 

I I 2 PITT 



7^2 

c 



$ ( J] - 2vrsign(fi) 



v(Q)dQ. 



r — 1 



s + 1 



s = l,...,r (86) 



1 2/ 1 < rn. (87) 



v(Q)dQ 



(88) 



Let us recall that the variance of the discrete white noise is a 2 = Additionally, let us define the function 



rN n 



T c J T c Efy 



v(tt). 



(89) 



By substituting (1191 ) in (1871) and (1881) . using definition (1891 ), and to — |f we obtain the fixed point equation 
(l26l) and the limit (1251) . respectively. This concludes the proof of Corollary El 



Appendix IV 
Proof of Theorem [2] 

The proof of Theorem [2] follows along the line of the proof of Theorem [Q In this case A.^ j7 .(£l, r) 
c — ^* / u i e ^-j anc i me ma t r ix r) is independent of r. Specifically, 



Q(^r) = — 



e{Q)e H (Q). 
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Then, applying the same approach as in Theorem Q] Lemma [5] and Lemma [6] yield 



lim SINR A 



W\ 2 r 
ItxT} 



with 



T c 



e(Q)e H (n) 



+oo 



e H (Q)T{a 2 ,Q)e{Q)dQ 



AdF| j4 | 2 (A) 



1 + 2w f-TT T? 



a 2 I r + (3 



U(Q)D'(tt)U H {tt)\dF lAl 2{\) 



1 + 2vr J-tt T} 



e H {W)T{(T 2 ,n')e(n')dn> 



e H (Q')r(a 2 ,Q')e(Q')dQ' 
(90) 



with C7(fi) defined in (l6Tl) and ZD'(O) diagonal matrix with all zero elements except the ([^J + l) 



th 



Then, it is apparent that the 



element, corresponding to the eigenvector e(Q) and equal to ^ 
solution of the fixed point matrix equation (l90l) is a matrix with the basis of eigenvectors U(Q), and (190 
reduces to the equation corresponding to the ( + l) th element v(Q) of T(fi) = {/^(fi) Y(f2)L/"(n) 

AdF| Ap (A) 



t;; 1 (fi)= ( r 2 + /3- 



(91) 



The other components of the diagonal matrix Y(f2) are simply given by t'7 1 (fi) = a 2 , s = l,...,r and 
s ^ ([EzlJ + l) . The identity e H '(fi)T(fi)e(fi) = yields 



lim SINR fc 



|qfc| 2 y 

2ttT 2 



u(fi)dfi. 



(92) 



The convergence (1921) in probability or in the first mean can be proven as in Theorem [T] By substituting (fT9~l) 
in (19TI ) and (|92l , using definition (|89l ), and c<j = |r, we obtain the fixed point equation (|36l) and the limit 
(1351) . respectively. 

This concludes the proof of Theorem [2l 



Appendix V 
Proof of Proposition [2] 

Proposition [2] follows immediately from Corollary [3l In fact, from (1331) it is apparent that the multiuser ef- 
ficiency of a system with load j3 and sine pulses having roll-off equal to a is equal to the multiuser efficiency 
of a system with load £ and sine pulses having zero roll-off. Thanks to the fundamental relations between 
multiuser efficiency and capacity ||39ll we obtain (|41"I) . Since the spectral efficiency is obtained as the ratio 
— , it is apparent from (1411) that it is constant as (3 — > oo for any finite bandwidth a. 



a 
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Appendix VI 
Proof of Corollary 

Let MMSE bfc [ m ](p) be the achievable MMSE by an estimator of the symbol bk[m] transmitted by user k 
in the m-th symbol interval when the transmitted signal b in ([141) is Gaussian and the signal to noise ratio 
is p = cr~ 2 . Furthermore, let SINR;, fc [ m ](p) be the SINR at the output of the same MMSE estimator for the 
transmitted symbol b k [m] . Then, 

MMSE «M = 1 + SINR, [m| W <93> 
Additionally, let I(b;y, p) be the mutual information in nats between the input b and the output y. From 

Theorem 2 in [[391 the following relation holds 

^-I(b ]y ,p) = E{\\Hb-Hb\\ 2 } (94) 
dp 

being b the conditional mean estimate. We recall here that for Gaussian signals conditional mean estimate 
and MMSE estimate coincide (see e.g., (351) and 

E{\\Hb - Hb\\ 2 } = tr (a 2 H H (HH H + a 2 iy l H) 



y- SINR bfcW (p) 
^p(l + SINR bfcW (p))- 



m,k 

For K,N,m — > oo with ^ — > {3, SINR bfe r m i(p) converges with probability one to deterministic values. 



More specifically, SINR fefcM (p) = ^L_^ 



being r/ 7 the multiuser efficiency corresponding to 7 



N ^7 

as in Corollary [2] or Theorem [2] 

Then, the total capacity per chip constrained to a given chip pulse waveform is given by 



ln2 7 J Q l + \tr) t 



This concludes the proof of Corollary @] 



Appendix VII 
Proof of Theorem [3] 

In this section Ojy denotes an A r x M matrix of zeros. Shortly, 0^ and denote A^-dimensional vectors 

(r) 

of zeros and ones, respectively. Additionally, we introduce the notation F N = <g> J r , with .F^ already 
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A basic property of the Fourier eigenbasis functions is stated in the following. 



41 



Property 1 Let E L denote an L x L matrix of ones and by JF^ (it), with u < (L — l)Nr, the LNr x LN 
matrix with structure 





^ r\ur 
^ur 


r\Nr 
™ur 




l)N 




nur 
u Nr 


F (r) 
1 N 


U Nr 


l)N 




nur 

\ U (L-l)N-u 


nNr 

U (L-l)A^u 


U (L- 


l)N 
1)N 



Then, 



Vl 



E L ®FP)F { [\u) = £(u,r,N) 



(98) 



where £(u,r, N) is a matrix with structure 



£{u,r,N) 



( 







LNr 



n 



rN 



7L 



rN 



f.(L~l)Nr-ur 
U (L-l)7Vr-ur 



and TZrN zs an r ^ x N block diagonal matrix with t-th block (71% 



rN) 



Let us consider the virtual spreading sequence of user k for symbol m in the time interval [—MT S , MT S 
with M > m integer: 



^ — a k[^(M+m)Nr^ki=k i ^ {M-m~l)Nr\- 



(m) r.T 



Property [Hand decomposition (|72|) yield 

1 



E 2M+ i ® F% ] )^ ] = e 2M +i ® (n T r k N A^ r ( Tk )s k 



V2M+ 1 



{E 2M+1 ®F^)i)t ] 



V2M+T 



(99) 



with 1)1'"' = a k [0f M+m)Nr , {F < pTZ T r k N A^ r (T k )s k ) T ,0j M _ m)Nr } T . Let us observe that the position of the 
nonzero elements does not depends anymore on f. Since the random entries of the vector s k are rotationally 
invariant, the rotation matrix TV rN can be absorbed in the random vector !> k without change of the statistics 
and the vector \] k can be rewritten as 

= a k[0(M+m)Nr, 0(M-m)Nr] T 



November 28, 2009 



Submitted manuscript to IEEE Transactions on Information Theory 42 

where hf^ = F$ A^ iT .(r fc ) is the virtual spreading for the m-th transmitted symbol of user k delayed by 
t G [0,T C ). From the previous considerations it follows that the random matrix T = 7t7i H is unitarily 
equivalent to an infinite block diagonal matrix with blocks T^ m ' = H^ m 'H^ m ' H of dimension rN x rN 
being the matrix with k-th column equal to fijj. , i.e., the transfer matrix of a symbol synchronous but 
chip asynchronous system with time delay {ri, . . . , tr}- Then, asymptotically for K,N — > +00 the eigen- 
value distribution of the matrix T^ m ' equals the eigenvalue distribution of the matrix T. The equivalence 
of the systems in terms of SINR at the output of a linear MMSE detector and in terms of capacity follows 
also from the fact that the SINR is invariant to any unitary transform of the system transfer matrix as already 
observed in the proof of Theorem [IJ 
This concludes the proof of Theorem [3] 
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